Dark Side of Higgs Diphoton Decays and Muon g — 2 
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We propose that the LHC hints for a Higgs diphoton excess and the muon g — 2 (gr M — 2) discrepancy 
between theory and experiment may be related by vector-like "leptons" charged under both U(1)y 
hypercharge and a "dark" U(l)d- Quantum loops of such leptons can enhance the Higgs diphoton 
rate and also generically lead to U{1)y-U{1)cl kinetic mixing. The induced coupling of a light U(l)d 
gauge boson Zd to electric charge can naturally explain the measured g M — 2. We update Zd mass 
and coupling constraints based on comparison of the electron g — 2 experiment and theory, and find 
that explaining <? M — 2 while satisfying other constraints requires Zd to have a mass ~ 20— 100 MeV. 
We predict new Higgs decay channels "fZd and ZdZd, with rates below the diphoton mode but 
potentially observable. The boosted Zd — > e + e~ in these decays would mimic a promptly converted 
photon and could provide a fraction of the apparent diphoton excess. More statistics or a closer 
inspection of extant data may reveal such events. 



The discovery of a new Higgs-like state, at a mass of 
about 125 GeV, by the ATLAS PJ and CMS @ col- 
laborations at the Large Hadron Collider (LHC) repre- 
sents a historical breakthrough in particle physics which 
is likely to provide a major step toward understanding 
electroweak symmetry breaking (EWSB). It remains to 
be seen whether this new state is the long-sought Stan- 
dard Model (SM) Higgs or some variant of it. While 
addressing this important question requires more data, 
the current results by both ATLAS and CMS hint, at a 
roughly 2a level 0, 13] , that the diphoton branching frac- 
tion of the new state, which we will henceforth refer to 
as the Higgs boson H, seems to be a factor of ~ 1.5 — 2 
larger than the SM prediction [!, H|. If that diphoton 
excess is confirmed with more statistics it would be an 
important clue for physics beyond the SM. 

Besides the recent hints from the Higgs data, the cur- 
rent 3.6cr discrepancy between the SM prediction and the 
measured value of the muon g — 2 denoted by 0^ — 2, 

is another potential clue pointing to new physics [7|. In 
this work, we propose that the Higgs diphoton excess and 
the g M — 2 discrepancy can be naturally related through 
the introduction of heavy new vector-like leptons charged 
under both U(l)y as well as a new U(l)d gauge symme- 
try with a relatively light O(20 - 100 MeV) Z d boson [|[. 
(The lower bound of 20 MeV will be discussed later.) We 
will refer to the new quantum number as "dark charge" 
Qd since Z& does not have direct couplings to the or- 
dinary "visible" SM sector. However, our study needs 
not assume a specific connection with dark matter (DM) 
physics, although it is a possibility. This issue will be 
briefly addressed. 

To set the stage for our discussion, we first address 
the main problem with trying to explain an excessive 
Higgs diphoton, H — > 77, decay rate. The SM predic- 
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tion arises from destructive interference of W ± loops (the 
dominant contribution) and a smaller (relatively nega- 
tive) top quark loop. Adding new heavy charged chiral 
fermions leads to additional negative loop contributions 
which further reduce the diphoton Higgs decay ampli- 
tude. (We do not consider cases where many fermions 
are added which change the overall sign and magnitude 
of the i?77 amplitude.) However, if the new fermions are 
vector-like doublets and singlets, the existence of heavy 
gauge invariant masses combined with mixing induced by 
the Higgs-singlet-doublet Yukawa couplings can change 
the sign of the new fermion loop contribution and actu- 
ally enhance the Higgs diphoton branching ratio. Vari- 
ants of that possibility have been suggested by a number 
of authors |9l-friT] who have discussed such scenarios in 
detail, including experimental constraints on properties 
of new vector-like fermions 18]. 

Here, we assume the above solution to the diphoton 
excess as our starting point, but endow the new fermions 
with an additional U(l)d gauge symmetry with a light 
gauge boson, Zd of mass mz d — 20 — 100 MeV. To avoid 
changing the Higgs production rate through gluon fu- 
sion, thereby affecting rates for other final states, these 
fermions should not carry SU(3)c color quantum num- 
bers. Hence, we focus on new "charged leptons." 

One-loop diagrams involving the new vector-like lep- 
tons can also induce, via kinetic mixing, a suppressed 
coupling of Z d to ordinarily charged particles. Such a 
gauge boson has recently been invoked in generic discus- 
sions of DM particles and their potential phenomenology 
[IjJIItJ. As we shall see, the Z d provides a natural viable 
solution to the — 2 discrepancy for a narrow range of 
mz d and leads to predictions for new Higgs decay modes, 
H — > jZd and ZdZd at the LHC, which may occur at 
observable rates or if not seen, used to constrain such 
models. 

Vector-like heavy fermions can be found in various new 
physics models including composite Higgs }2l| and super- 
symmetric U(l)' models [22j motivated to address the \i- 
problem [23| . To avoid problems with new stable charged 
particles, we consider charges which are multiples of the 
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FIG. 1: New fermion (F) loop contribution to H — > 77 decay 
(left) and 7 — Z<j kinetic mixing (right). 



SM particle charges. This would typically imply that 
the new particles carry hypercharge (for the sake of min- 
imality, we do not consider SU(2) triplet fermions as they 
would require additional Higgs content). 

Let us denote the ratio of the enhanced rate for H — > 
77 compared to that in the SM by 



R 
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Using the results of Ref. [1J|, the contribution of a new 
fermion F of electric charge Q and mass mp (see Fig. [T]) 
to the above ratio is given by 



(4/3)Q 2 dlogm F 



-4; 



SM 
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where Asm — —6.5 stems from the SM amplitude for the 
decay, v = V2(H) ~ 246 GeV, and m H ~ 125 GeV is 
the Higgs mass. Eq. ((5J) implies that for i? 77 > 1 we need 
9 log nip jd log u < 0, that is, the contribution of EWSB 
(v 7^ 0) to the mass of F must be negative. 

We next describe a simple scenario which explains the 
diphoton excess. We extend one of the examples pro- 
posed in Ref. [3] to include dark U(l)d interactions. The 
extended model contains vector-like "leptons", ip and x 
with SU(3) C x SU(2) L x U(1) Y x U(l) d charge assign- 
ments 

(^,V c )~(l,2) {±i±1} ; (x,X c )~(M) {t i, t i } , ( 3 ) 

where upper (lower) signs are for ifi (ifi c ), etc. Here, elec- 
tric charge Q is related to hypercharge Y and the third 
component of isospin T3 by Q — T3+Y. The mass matrix 
is obtained from the following Lagrangian: 



which is valid for y ~ y 

2 yy c v 2 



with 



A- 



and A, 
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The enhancement i? 77 ~ 1.5 can be achieved for yy c ~ 1 
and mi ~ m#, assuming b? m -C A 2 . 

Without further assumptions, the model in Eqs. ^ 
and (|1]) will lead to stable charged particles since the 
lightest new vector-like lepton L\ will be char ged . To 
avoid this unphysical situation, following Ref. we 
mention two approaches: 

(I) One possibility is to allow very small mass mixing 
between the new charged heavy leptons Li (i — 1,2), 
and ordinary SM charged leptons £j (j = e, /1, t) which 
gives rise to Li — > £jZd decays. For us, this could be 
simply achieved if we introduce a "dark" Higgs field, (fid, 
with Qd = — 1 which allows small Yukawa interactions 
(1) Gij<fi>dL\lj where Gij is a general Yukawa coupling matrix 
which connects the new heavy charged leptons with or- 
dinary charged SM leptons. The left-right and right-left 
Gij couplings will in general be different. The (fid is well 
motivated since it can also provide a mechanism for spon- 
taneous symmetry breaking and lead to mz d < 100 MeV 
if {(fi d ) ~ 100 MeV K ~ 9d <^». 

The Yukawa mixing interaction will induce potentially 
dangerous flavor changing weak neutral current inter- 
actions for the Z, H and Zd bosons. Here, we focus 
on the Zd couplings. They are important because the 
small mz d leads to interesting enhancements and are spe- 
cific to the model we are considering. The induced non- 
diagonal interaction Gij(mz d /mL i )Z^£jj fl Li appears to 
be highly suppressed by the mz d /mi Ji < 10~ 3 factor. 
However, for the longitudinal component (or Gold- 
stone mode s d ), that factor is cancelled 24 1 and one 



finds the coupling ~ GijS d £jLi as required by the Gold- 
stone boson equivalence theorem with different Gij for 
left and right handed Li. To avoid generating large chi- 
ral changing loop effects in quantities such as the electron 
and muon anomalous magnetic moments, lepton number 
flavor violating amplitudes for [i — ¥ ey, r — ¥ ^7, light 
lepton loop induced masses, etc, some of the \Gij\ (par- 
ticularly those involving e or jj) have to be quite small 
< 10~ 3 . However, even with such small couplings, the 
decay rates T(Li — > tjZd) ~ IGjjpmi^ are likely to pro- 



-C m = m^ifiip + m x xx c +yHifix+y c ^ r V^x' : +fi-c (4) vide very prompt signals, Li -» £jZ d , at the LHC, where 

the light Zd, which subsequently decays into e + e _ , can 



After EWSB, we get two Dirac fermions L\ and L2 with 
electric charge \Q\ = 1 and masses m\ and m%, where 
mi + 7712 = Ttity + rn x , 7772 > m\ and one neutral Dirac 
fermion N of mass mj, with mi < < m 2 . All the 
new heavy leptons have pure vector couplings to gauge 
bosons; hence, no gauge anomalies are present. 

One can show that the above field content results in 
an enhancement of the diphoton rate, which to a good 
approximation is given by [14| 



R 



77 



1 + 0.1 



1 + v/Af + 



(5) 



mimic a converted high energy photon. (One expects 
LiLi pairs at colliders actually giving rise to di-£ jZd 
events.) We note that by assuming small mixing param- 
eters (< 10~ 3 ), we can avoid conflict with precision mea- 
surements of e-fjj-T universality such as those discussed 
in Ref. [25j]. Of course, one might use lepton mixing ef- 
fects to make further predictions or to accommodate all 
or part of the g M — 2 discrepancy. The constraints on this 
model and its phenomenology are potentially rich and 
interesting, but beyond the scope of this paper. Here, we 
only suggest it as a means to avoid stable heavy charged 
leptons. 
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(II) Another possibility is to avoid mixing with the SM 
leptons and instead add the fields (n, n c ) ~ (l,l)r 0= pu 
(carrying only Qd ^ 0), with the new expanded neutral 
fermion mass terms 

m n nn c + y n H^n + y c n H^) c n c + h.c., (7) 

which together with Eq. (j4j result in two neutral Dirac 
fcrmions, n\ and ri2, with masses m„ 2 > m ni and po- 
tentially m ni < mi. Given the mixing terms in Eq. ([TJ, 
the new charged L\ particle would decay into (possibly 
virtual) W^- and n%. 

It may be tempting to think of the lightest new neu- 
tral state, n\ as a stable relic DM candidate. However, 
without further assumptions, this turns out to be not 
phenomenologically viable. We will address this ques- 
tion and discuss potentially viable DM alternatives in 
the appendix. 

We now turn to the g M — 2 problem which has persisted 
over the last several years. The discrepancy between the 
measured value and the SM prediction is about 3.6c [5L 

Aa M = a)f p - a™ = 287(80) x 1(T 11 , (8) 

where a M = (g^ — 2)/2. A simple explanation [26| of 
this difference postulates a new hidden U(l)d symmetry 
which kinetically mixes with £7(1) y by 

where e parametrizes the mixing, 9yy is the weak mixing 
angle, and X^ v = d^X v - d v X^ X = B,Z d , is a {7(1) 
field strength tensor. Upon kinetic diagonalization, one 
finds that the massive dark boson Zd obtains an induced 
coupling eeZjJ^™, where J^ m is the the electromagnetic 
current (27J- This light boson is a target of active and 
planned searches at JLAB and MAMI in Mainz [29M35j . 
Early results from those experiments are illustrated as 
constraints in Fig. [5] 

One can show that the 1-loop contribution of the Zd 
to a M is given by [26|, |36[ 

a** = -^e 2 F v (m Zd /mJ (10) 

withFv(ar) = dz[2z(l - z) 2 ]/[(l - z) 2 + x 2 z] ; F v (0) = 
1. 

In Fig. [21 we give the current exclusion bounds on e 2 
(adopted from Refs. [28|, |3l|). There, we have updated 
the bounds coming from the recently improved electron 
anomalous magnetic moment comparison between exper- 
iment H3 and SM theory [Mill: 

Aa e = < xp - a s e M = -1.06(0.82) x 10~ 12 . (11) 

Because of the small momentum transfer in Rydberg 
measurements Q 2 <C m| d , the effect of a light Zd on the 
determination of a in Ref. [38| is expected to be negligi- 
ble for the m Zd mass range considered. That constraint 




m zd [MeV] 

FIG. 2: Exclusion region in rnz,— e 2 space, updated from 
Refs. 0, (3l[ to include recent [33-[3§] a e results ("new a e "); 
for comparison we also present the old bound ("old a e "). 

implies at the 3tr level af d < 1.4 x 10~ 12 (with af d ob- 
tained from Eq. (fTOf with m M — > m e ) which rules out a 
significant region which would otherwise provide a viable 
Aa^ solution (i.e. m Zd < 20 MeV, e 2 - 0{2 - 5 x 10~ 6 ). 
Hence, one finds that the discrepancy in Eq. ([5]) can be 
explained for 

20 MeV < m Zd < 100 MeV (12) 

and 

2 x 10~ 6 < e 2 < 10- 5 (13) 

without conflict with current experimental bounds on Zd- 
We note that the exclusion region due to a e is somewhat 
enhanced because of the sign of Eq. (fTTj) which is opposite 
to that expected from Zd- 

Our new constraint rules out the (previously allowed) 
10 MeV < m Zd < 20 MeV region of the "a M " band in 
Fig. [2j That explicit part of the band is the focus of a 
proposed direct Zd search at VEPP-3 [4(J; however, the 
experiment will also explore smaller e 2 . 

In a simple U(1)y x U(l)d framework, e is an arbitrary 
renormalized parameter set by experiment. Normally, we 
expect e ~ egd/8n 2 which for gd — e gives e ~ 10~ 3 . 
That expectation would be natural, if either of the U(l) 
symmetries at low energy descend from a non-abelian 
group in the ultraviolet sector, such that the high energy 
(bare) value of e is zero. In that case, finite kinetic mix- 
ing of the type in Eq. ^ can be naturally induced by 
loops of fermions which are charged under both U(1)y 
and U(l)d [4l|. The typical value of the kinetic mixing 
at low energies can then be estimated from a 1-loop dia- 
gram (see Fig. [T]) , which is roughly in the range indicated 
by Eq. (|13|) . In fact, given an appropriate assignment of 
fermion charges, one can calculate a finite 1-loop result 
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FIG. 3: New fermion (F) loop contribution to H — > jZd. 
decay (left) and H —¥ ZdZd decay (right). Both charged and 
neutral fermions should be included in the loop on the right. 



[41j . For example, a 2-generation extension of the refer- 
ence model in Eq. ([3]) with opposite sign dark charges Qd 
results in a finite and computable value 



67T 2 



log 



\m 2 m 3 ) 



(14) 



where 7713,4 correspond to the masses from the second 
generation of charged vector-like leptons. We see that 
for typical values of parameters, e in Eq. (fl4"|) has the 
size needed to address the g^ — 2 anomaly. For example, 
if the logarithm involving the masses is order unity, for 
Qd ~ 1 and gj ~ e we find e ~ 1.5 x 10~ 3 . So, the model 
introduced in Eq. (|3j to explain a H — > 77 excess can 
also accommodate the g^ — 2 discrepancy. 

Our model leads to the interesting prediction that the 
Higgs has new decay channels yZd and ZdZd (Fig. [3J 
with rates somewhat smaller than that of the 77 mode. 
To see this, note that Eqs. (fl"3|) and (fT4|) imply we need 
gd ~ e to explain the measured g^ — 2. Hence, we expect 
that the new decay modes in Fig. [3J will have a similar 
amplitude as the extra contribution to 77 (Fig. [I]). 

To connect the Higgs diphoton excess and g^ — 2, as 
proposed here, it is sufficient to have a single fermion 
which carries both U(l)y and U{\)d charges. This min- 
imal setup can effectively emerge in our reference model 
in Eq. (J3J if there is a modest hierarchy of masses and 
the new Higgs decay amplitudes (Figs.Q]and[3]) are dom- 
inated by the lightest charged state. Under such a sim- 
plifying assumption, a rough estimate for the rate of the 
H — > yZd compared to the observed rate of H — > 77 (not 
the SM expectation) can be given in terms of i? 77 by 



T(H -» yZ d ) 
Y(H -> 77) 



2 1 



1 



(g± 

\ e 



~ ■ (15) 



where the factor of 2 accounts for the nonidentical final 
state particles, and the first set of parentheses factors 
out the new lepton contribution to the H ~ > 77 decay. 
Similarly, the rate of the H — » ZdZd compared to H — > 
77 is 



r z d z d 



_ T(H -> Z d Zd) 
T(H -+ yy) 



ay 



( 16 ) 

Based on our preceding discussion, let us take g& = e 
as a typical value for our scenario. This would im- 
ply r lZd « 0.07 - 0.17 and r ZdZd w 0.03 - 0.09 for 



0.20 



0.15 



0.10 



0.05 



().()( 



.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Rile. 



FIG. 4: r 1 z d (rz d z d ) bounded by blue solid (red dashed) 
curves versus gd/e, for _R 77 = 1.5 (bottom) — 2.0 (top). 



i? 77 = 1.5 — 2.0 (see Fig. @|. We see that the rates for 
these new channels are expected to be well below that of 
H ^ 77 but potentially within the reach of the LHC ex- 
periments. Note that for mz d in the range of Eq. ( TT21) . the 



Zd will mainly decay promptly into e + e~ [24|, i.e. within 
the beam pipe region. However, since mz d *C tuh, the 
decay products will be boosted and highly collimated. 
In general, we expect that the Zd in the final state will 
mimic a promptly converted photon (7 — > e + e~) but 
with a small nonzero mass and production vertex near 
the beam rather than in the tracker. Such properties 
could be used as a signal for these new decays. 

For the m Zd = 20 — 100 MeV range, the opening an- 
gle of order 0.002 is well below the 0.02 roughly required 
for separating the e + e~; so, except for the effect of the 
magnetic field in the detector, the e + e~ pair would be 
indistinguishable from a photon (see Ref. [42| for some re- 
lated discussions). However, the magnetic field will lead 
to separated tracks for e~ and e + in the tracker. Un- 
der the assumption of a converted photon, these tracks 
would normally be fitted assuming zero invariant mass. 
For our purposes, it would be useful to modify the fitting 
program to allow masses 20 — 100 MeV. A careful exam- 
ination of the available diphoton data, which is outside 
the scope of this work and more appropriate for exper- 
imental scrutiny, could reveal or constrain the presence 
of H — > yZd or even H — > ZdZd events. 

The above rough estimates of the relationships between 
77, yZd and ZdZd decay rates can change if the under- 
lying model deviates from the simplifying assumptions 
used to derive Eqs. (fT3|l and (fTB|) . For example, in a 2- 
generation extension of the model, the charged leptons 
carrying opposite dark charges are likely to cancel par- 
tially and somewhat reduce the rate of H —> yZd (Fig.[3j. 
In this circumstance, our estimate of the H — > yZd 
should be viewed as an approximate upper bound. In 
the case of the H — > ZdZd amplitude (Fig. [3J, if the 
model is enlarged to include the interactions in Eq. (J7J, 
the resulting neutral leptons could either increase or re- 
duce the rate for H — > ZdZd- In addition, a light Zd may 
be only a small part of the Aa M discrepancy, consistent 
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with a much smaller e 2 below the "a M explained" band 
in Fig. [3J In that case, the jZj and rates could be 

much reduced. 

Overall, our estimates for the new H decay rates are 
meant to be suggestive and to stimulate experimental 
searches for those decays. Definitive predictions require 
more detailed studies using specific models and param- 
eters. We note that the search for H — > ^Zd and 
H — > ZdZd with Zd resembling a promptly converted 
photon in the LHC experiments would be largely com- 
plementary to light "dark boson" searches such as those 
at JLAB and in Mainz (29l - [35| . Those programs will not 
only directly probe the a M explained region in Fig. [2] for 
Zd boson, but will also explore significant parts of pa- 
rameter space outside that band. 

"Dark" decay modes of the Higgs may also arise 
through other mechanisms 0, ED, 0] • For example, it 
is possible to have a H — > ZdZd mode from Higgs (H) 
and dark Higgs (fa) mixing [43|, and in the presence of 
mass mixing between Z and Zd, as studied in Ref. [24j . 
a new Higgs decay mode H — > ZZd would also be possi- 
ble. The latter channel could mimic Z7 with a promptly 
converted photon if Zd is sufficiently light. However, the 
predictions in those cases are more arbitrary; whereas the 
connection between the Higgs diphoton rate and g^ — 2 
in our model allows us to make more quantitative esti- 
mates for the dark decay rates of the Higgs. We also note 
that our loop induced "fZd and ZdZd H decays involve 
primarily transverse Zd bosons while the ZdZd and ZZd 
decays in Refs. [24j,|43[ are dominated by longitudinal Zd- 

In this paper, we have discussed a possible link between 
the reported excess of the Higgs to diphoton decay at the 
LHC experiments and g^ — 2 via heavy new vector-like 
leptons and a light dark gauge boson. A gauge boson 
of mass ~ 20 — 100 MeV with small induced coupling 
to the SM particles is well motivated as a rather simple 
explanation of the 3.6<r deviation of g^ — 2 from the SM. 
The required coupling of the Zd to the SM fermions is 
naturally obtained when it arises from loops of charged 
extra fermions that couple to both the SM U(l)y and 
a dark sector U{l)d with similar size couplings. This 
scenario yields an additional contribution to H — > 77 
through a loop of the charged extra fermions, which is 
consistent with the recent 2a level deviation at the LHC 
experiments. The Higgs boson can also decay into jZd 
and ZdZd, with the light Zd bosons looking like promptly 
converted photons in the ATLAS and CMS detectors. 
Such a connection implies a few definite predictions in 
high energy experiments at the LHC and complementary 
low energy Zd searches at JLAB and in Mainz. 
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Note Added: After this paper was posted and sub- 
mitted for publication, a preprint (47j appeared which 



reached similar conclusions regarding the u pda ted exclu- 
sion region due to new a e results in Refs. |38l |39| . The 
authors of Ref. [13] also gave a detailed analysis of the 
effect of a Zd on a explicitly showing that for the mz d 
range considered, it has a negligible effect. 



Appendix A: Relation to Dark Matter 

As mentioned before, the lightest neutral Dirac lepton 
in our simple framework, as presented above, is not a 
good DM candidate. Starting from our basic model in 
Eq. (J3J, let us examine solutions (I) and (II), discussed 
earlier, for avoiding a stable charged state L\. In case (I), 
a small degree of mass mixing with the SM leptons would 
not alter the spectrum of the model significantly, and we 
still expect the neutral state N to be more massive than 
Li and thus allow the decay N — > L\ + (virtual) W. 
Hence, N cannot be a long-lived DM candidate. 

In case (II), the neutral particle ri\ is the lightest new 
state and could be stable. However, it carries dark charge 
and can hence elastically scatter from protons, through 
light Zd exchanges, at too high a rate. For the param- 
eter space of interest here, we would typically expect 
the scattering cross section of n\ from a nucleon to be 
~ 10~ 33 cm 2 (see, for example, the last paper in Ref. [l9j]) 
which, for masses ~ 100 GeV, is ruled out by about 11 
orders of magnitude from direct DM search constraints 
45]. 

One could arrange for the early universe annihilation 
cross section of n\ to be large, so that its relic density 
is very suppressed and it is not a primary component 
of DM. Alternatively, one could instead take a positive 
approach and extend our framework to allow for the pres- 
ence of a DM candidate. For example, let us postulate 
the U(l)d breaking scalar fa, with (fa) ~ 100 MeV, 
which was invoked earlier in our discussion of case (I). 
We also add a singlet vector-like lepton £ without any 
gauge charges and assume that all of the new fermions 
are odd under a Z 2 parity to forbid DM decay. With £ in 
the spectrum, we can write down an interaction A0dn£, 
with A ~ 1. If m ni > TOf, the above interaction would 
lead to tt-1 — > fa £ which will be prompt, and all m parti- 
cles will decay into £ and fa. The scalar <pd will eventu- 
ally decay into SM states. However, £ would be a viable, 
stable DM candidate. Note that due to the very small n £ 
mixing induced by (fa) ^ 0, the interactions of £ with 
Zd could be quite suppressed, and hence one would avoid 
the stringent bounds from direct detection. 

The relic density of £ is set by £ £ — > fa fa, through t- 
channel exchange of neutral states, which for weak scale 
masses could be of the correct thermal relic size. The 
n£ mixing is of order (fa)/m n ~ 10~ 3 , assuming ~ 
100 GeV, typical of the particles considered here. Hence, 
the direct detection cross section via Zd exchange with 
protons could be suppressed by ~ 10 -12 , which is just 
below the current sensitivities [451 ] . 

An alternative possibility entails adding a lepton num- 
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ber violating interaction to our model which splits the 
neutral Dirac lepton n\ and its antiparticlc partner into 
two nondegenerate Majorana states. Such states have 
zero dark charge and will not scatter elastically (at least 
at leading order) off ordinary matter via Zd exchange 



i46j. There will be off-diagonal Zd couplings between 
the two Majorana states which allow inelastic scattering, 
but for the lighter state, that can be kinematically sup- 
pressed. A detailed evaluation of this scenario is beyond 
the scope of our study. 
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